Mon. Not. R. Astron. Soc. OOP. [THTSl (2010') 



Printed 4 November 2010 



(MN I*TeX style file v2.2) 



The Rising Star- Formation Histories of Distant Galaxies 
and Implications for Gas Accretion with Time 



o 

> 

O 



Casey Papovich^'^*, Steven L. Finkelstein^'^, Henry C. Ferguson^, 
Jennifer M. Lotz^'^^f , and Mauro Giavalisco^ 

^Department of Physics and Astronomy, Texas A&M University, College Station, TX 775^5-4242, USA 
^George P. and Cynthia Woods Mitchell Institute for Fundamental Physics and Astronomy 
^ Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 
'^National Optical Astronomy Observatory, 950 N. Cherry Avenue, Tucson, AZ 85719, USA 
^Department of Astronomy, University of Massachusetts, Amherst, MA 01003, USA 



cn 



Accepted 2010 November 1. Received 2010 October 18; in original form 2010 July 12 



o 
u 

o 



(N 
> 

^, 

i> 

o 
o 



X 



ABSTRACT 

Distant galaxies show correlations between their current star-formation rates (SFRs) 
and stellar masses, implying that their star-formation histories (SFHs) are highly 
similar. Moreover, observations show that the UV luminosities and stellar masses 
grow from z=8 to z=3, implying that the SFRs increase with time. We compare the 
cosmologically averaged evolution in galaxies at 3 < z < 8 at constant comoving 
number density, n = 2 x 10~^ Mpc~'^. This allows us to study the evolution of stellar 
mass and star formation in the galaxy predecessors and descendants in ways not 
possible using galaxies selected at constant stellar mass or SFR, quantities that evolve 
strongly in time. We show that the cosmologically averaged SFRs of these galaxies 
increase smoothly from z = 8 to 3 as ^(i) '^ i" with a = 1.7 ± 0.2. This conflicts 
with assumptions that the SFR is either constant or declines exponentially in time. 
Furthermore, we show that the stellar mass growth in these galaxies is consistent with 
this derived SFH. This provides evidence that the slope of the high-mass end of the 
IMF is approximately Salpeter unless the duty cycle of star formation is much less than 
unity. We argue that these relations follow from gas accretion (either through accretion 
or delivered by mergers) coupled with galaxy disk growth under the assumption that 
the SFR depends on the local gas surface density. This predicts that gas fractions 
decrease from z = 8 to 3 on average as /gas ~ (1 + z)^-'^ for galaxies with this number 
density. The implied galaxy gas accretion rates at z > 4 are as fast and may even 
exceed the SFR: this is the "gas accretion epoch". At z < 4 the SFR overtakes the 
implied gas accretion rate, indicating a period where galaxies consume gas faster than 
it is acquired. At z < 3, galaxies with this number density depart from these relations 
implying that star formation and gas accretion are slowed at later times. 

Key words: cosmology: observations - galaxies: evolution - galaxies: formation - 
galaxies: fundamental parameters - galaxies: high-redshift - galaxies: stellar content 



1 INTRODUCTION 

The process by which galaxies acquire their baryonic gas 
and convert it to stars is one of the main questions for the- 
ories of galaxy formation. However, it has been very diffi- 
cult to constrain directly the evolution of these processes 
primarily because it is challenging to link the descendants 
and progenitors of galaxies at different snapshots in redshift. 
Currently, we have only a weak observational understand- 
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ing of the star-formation histories (SFHs) of distant galax- 
ies, which is arguably one of the most important quantities 
to measure because it governs directly how rapidly galaxies 
build up their stellar content. 

Nevertheless, there are recent hints from the interpre- 
tation of observational data that distant galaxies experience 
SFHs where the SFR increases with time (decreasing red- 
shift). The relatively tight correlation and low scatter be- 
tween the star format ion rates (SFRs) and stellar masses in 
distant galaxies (e.g..lDaddi et al.ll2007l : iNoeske et"ai]|2007l : 
IStark et al1l2009l : lMagdis et al.ll2010l ') imply that the major- 
ity of galaxies at high redshift sustain high levels of star for- 
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mation for prolonged time-scales of ^^0.5-1 Gyr, with high 
"duty cycles" |j Therefore short-lived starbursts driven by 
major mergers likely play a minor role. If the SFRs of these 
galaxies have been constant or declining in time, then it im- 
plies that the progenitors of these galaxies would have had 
higher specific SF Rs (SFRs per unit stellar mass). However, 
IStark et al.l ([20091) observed a lack of evolution in the specific 
SFRs of UV-selected galaxies from z = 4 — 6, which excludes 
the possibility that galaxies at z = 4 have been forming stars 
at con stant or declining rates since 2; > 5. iFinkelstein et al.l 
l|2010l) noted that both the characteristic luminosity of the 
UV luminosity function, Luvi ^-nd the stellar mass of galax- 
ies at the characteristic luminosity increase with decreasing 
redshift from z ~ 8 to 3, which is most easily explainable if 
distant galaxies experienced SFRs that increa s e wit h time. 
Furthermore, recent work by iMaraston et al.l (|201G| ) argue 
that rising SFHs better reproduce the observed colours of 
galaxies at z ~ 2 compared t o models that decline with 
time (see also. lLee et al.ll2010bl ). 

Support for periods of extended star formation, such 
as that expected from a rising or constant SFH, also 
comes from kinematic measurements of many star-forming 
galax;ies at z ~ 1.5-4, which show evidence for extended 
rotating gase ous disks no t domi nated by spheroids or 
mergers (e.g.. iGenzel et al.1 ^2008: Da ddi et al.ll2009l. l2010al: 
Forstcr-Schrcibcr et al. 2009; Law ct al. 2009: ' Wright et all 
2009 ; Carilh et al. 2010; T acconi et al. 2010 ). The existence 



of gaseous disks likely requires that distant galaxies form 
stars at a quasi-steady pace, and not in short periods of en- 
hanced star formation s uch as merger-driven st a rbursts (see, 
Elmosrcon et al. 200»; iBournaud et all l2008l : iDekel et al.) 
2009 : Renzini 2009 ') . Comparing galaxy SFRs to their gas 



masses yields gas consumption time-scales of a few hundred 
Myr, typically much shorter than a Hubble time at the ob- 
served redshift. In order for the galaxies to maintain their 
high levels of star formation, it follows that they must re- 
ceive cold gas from the intergalactic medium (IGM) to re- 
plenish their su pplies of fresh gas to fuel the contin ued star 
formation ('e.g.. lErb|[2008l : |Prochaska fc Wolfell2009l l . 

SFRs that increase with time for z > 3 are an out- 
come from hydrodyn amic and semi-analytic simulations. 
IMaraston et al.l (|201G| ) favor SFRs that increase with time 
in part because these models reproduce the SFRs and stel- 
lar masses of galaxies from their s emi-analyt i c mode l. Sim- 
ilar conclusions are obtained by iLee et al.l ll2010bf) using 
SFHs from a semi-analytic model based on lSomerville et al.l 
1120011). Using a co smological hydrodynamic simulation, 
iFinlator et al.l ((20071) predict that galaxies at 2 ^ 6 expe- 
rience SFRs that increase with time, and they show that 
these SFHs reproduce the UV lu minosity functions and ob- 
serve d colours of distant galaxies (jFinlator et al.ll2006l . 120071 . 
[2OI0I). Interestingly, the predicted SFHs in these simulations 
are scale-invariant, varying in mass only by a scale factor, 
implying that measurements of the SFHs of distant galax- 
ies have strong constraining power on models. In the sim- 
ulations, these smoothly-rising SFRs are driven by the net 
gas accretion (gas inflow minus gas outflow) regulated by 



^ Here we take the duty cycle to be the ratio of the time duration 
of star-formation in galaxies relative to the total tim e interval 
probed by the galaxy survey, see, e.g., lLee et al.l 1120091) . 



the competition between the growth o f dark-matter haloes 
and decrease in c osmic density (e.g., iBouche et al.l |201G| : 
iDutton et al.ll201Gf ). Therefore, measuring the galaxy SFHs 
constrain both the growth of galaxy stellar mass, and the 
rate that gas accretion occurs. 

Here, we use observations of galaxies to measure em- 
pirically their SFHs from 2 = 8 to 3. We compare galaxies 
at constant (comoving) number density over a range of red- 
shift because this allows us to measure the evolution of the 
predecessors and descendants of galaxies in a relative mean- 
ingful way that is not possible using other selection methods. 
Our study is motivated in part by previous studies of galax- 
ies at constant number density to infer evolution (see, e.g., 
IConrov fc WechsleijIJOOgl : Ivan Dokkum et al.ll2010l ). Studies 
of galaxy properties at fixed mass or luminosity compare 
galaxies of very different number density as a function of red- 
shift, leading to potential biases (a variant of the well known 
Malmquist btas) when comp aring galaxie s selec ted at dif- 
ferent redshifts. For example, IStark et al.l (|2009l ) study the 
evolution of galaxies at constant rest-frame UV luminosity, 
finding negligible evolution in the average stellar masses and 
ages of galaxies over the redshift range 4 < 2 < 6. However, 
the number densities of galaxies at fixed luminosity change 
by factors >2 over this range (see figure[T]and discussion be- 
low) implying strong evolution. Comparing galaxy samples 
at constant number density has the beneflt that it allows us 
to follow the evolution of the progenitors and descendants 
of galaxies across different redshifts (see Appendix IX|| . 

This paper is organized as follows. In §[2] we discuss the 
properties of galaxies at constant number density at red- 
shifts 3 < 2 < 8. In § 131 we constrain empirically the SFHs 
of galaxies at constant number density. In § |4] we compare 
the joint evolution of stellar mass and SFR for galaxies at 
constant number density, finding them to be largely con- 
sistent. In § [5] we use the SFR and stellar mass evolution 
combined with the expected growth of galaxy disks to pre- 
dict the evolution of galaxy gas fractions at 3 < 2 < 8. 
We provide discussion in support of using galaxy samples 
at constant number density in Appendix [X] and we dis- 
cuss how models with SFRs that increase with time af- 
fect the stellar masses of distant galaxies in Appendix [B] 
Throughout this paper we use a cosmology with parame- 
ters. Ho = 70 km s"^ Mpc"\ r2A/,o=0.3, and Ao=0.7. All 
magnitudes quoted here are measured with resp ect to AB 
system, mAB = 31.4-2.5 log(/^/lnJy) (jOke fc Gunn.1983, ). 

2 THE PROPERTIES OF GALAXIES AT 
CONSTANT NUMBER DENSITY 

We compare galaxies at constant comoving number density, 
n = 2 X lO""* Mpc"^. For redshifts 3 < 2 < 8 this corre- 
sponds approximately to the space density of >L* galaxies, 
and at 2 ~ 0.1 galaxies with this number density have stel- 
lar masses >1.5 x 10^^ Mq, approximately twice that of the 
Milky Way Galaxy. Therefore, this number density corre- 
sponds to "typical" galaxies at these redshifts. 

Comparing galaxies at constant number density has sev- 
eral advantages. In the absence of mergers the comoving 
number density of galaxies is invariant with time: galaxies 
will grow in baryonic mass through gas accretion and grow in 
stellar mass through star formation, but their number den- 
sity will remain unchanged. In principle, identifying galaxies 
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at constant number density directly tracks the progenitors 
and descendants of the galaxies at different redshifts. Even 
a non-negligible number of mergers would have only a small 
effect on the progenitors and desc endants of galaxies selecte d 
at constant number density fe.g. Jvan Dokkum et al.ll201Cr ). 
In Appendix [A] we study how well selecting haloes at con- 
stant number density identifies their progenitors and descen - 
dants using simulated merger trees (ISpringel et al.ll2005a ). 
We found that selecting haloes at constant number density 
at different redshifts identifies 60-80% of the halo progen- 
itors and descendants over the entire redshift range from 
3.0 < z < 7.3, and that this completeness fraction remains 
high even when selecting at constant number density based 
on SFR instead of halo mass. In contrast, studies of galaxies 
selected constant mass or luminosity correspond to highly 
heterogeneous samples at different redshifts for the reason 
that these quantities evolve strongly with time. Therefore, 
selecting galaxies at constant number density is a relatively 
straightforward means to study the properties of galaxy pro- 
genitors and descendants at different redshifts in a way that 
is unfeasible using other galaxy selections. 

Here, we compare galaxies at constant number density 
based on their UV luminosity. Observations of gala:xy clus- 
tering show that gala xy UV luminosity traces the halo mass 
tight l y (e.g.,[G iavalisc o fc Dickinsorj|200ll : lAdelberger et al.l 
I2OO5I : iLee et all 12009 . 120091 ). Theoretically, the relation is 
unclear, and predictions span a wide range for the scat- 
ter between SFR (UV luminosity) and halo mass. While 
some models predict a ti g ht relation, with sca tter <0.1 dex 
(e.g., iDutton et aDbOld : iFinlator et"alll2010l'l. others pre- 
dict sizable s catter, 
iBertone et"!!.! 20071 



-2 



De Lucia etHI I2OO6I : 



~0.3 dex (e.| 

To avoid an over-reliance on theoreti- 
cal expectations, we make here the simple assumption that 
there is a one-to-one correspondence between the number 
density of galaxies at fixed halo mass and galaxies at fixed 
UV luminosity, which is supported by the observational evi- 
dence. However, in Appendix|X]we study how scatter in the 
relationship between UV luminosity and dark-matter halo 
mass impacts the completeness in the descendants and pro- 
genitors of galaxies selected at constant number density. To 
summarize, we find that increasing the scatter in this rela- 
tion reduces the completeness in the number of halo descen- 
dants identified at lower redshift. However, even in models 
where the scatter is highest, the completeness is >40% over 
3 < z < 7.3, implying that the method tracks the descen- 
dants with high fidelity. Furthermore, the majority of the 
descendants of galaxies that are "missed" at lower redshifts 
have halo masses within 0.3 dex of the selected galaxies. 
We find some empirical evidence that galaxies have aver- 
age SFHs that differ only by a scale factor related to total 
mass (§151), consisten t with expectations from simulations 
IFinlator et al.l (|2010|). Therefore, the descendants that are 
missed likely have similar SFHs. Even though the descen- 
dants of galaxies at selected constant number density are 
incomplete, this method allows us to study the evolution of 
galaxy populations over this redshift range. 



2.1 Selecting Samples 

Many studies have measured UV lu minosity functions of 
galax i es from z ^ 2 to z > 8 (e.g.. JBouwens et a.1.1 12004 . 
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Figure 1. The comoving spatial number density of galaxies 
brighter than given absolute magnitude, Muy, derived by inte- 
grating the UV luminosity functions at z = 3 to z = 8 using values 
in tablo[T]from Reddy & Steidel (2009) at z = 3.0, Bouwens et al. 
(2007) at z = 3.8 to 5.9, and Bouwens et al. (2010a) at z = 6.8 
and z = 8.0. We compared galaxies at constant number density 
by identifying at each redshift the absolute magnitude such that 
n(<M(jy) = 2 X 10"'* Mpc~'^ denoted by the dashed horizontal 
line. The arrows indicate the value of M^y that corresponds to 
this number density for each luminosity function. 



iMcLure et al.ll201Cl l. We integrated the UV luminosity func- 
tions to derive the number density of galaxies brighter than 
a given absolute magnitude at each redshift. 



n{<M) = 



(t>{M) dM 



(1) 



where (j>{M) is the luminosity function in units of number 
Mpc~^ mag~^. Figure [l] shows the integrated luminosity 
functions at each redshift. For the samples in this study 
we measured the evolving UV absolute magnitude, AIIjy{z), 
corresponding to the constant number density n(<i\fuv) = 
2 X 10"* Mpc~^ at each redshift. These values are given in 
table [T] and indicated in figure [T] 

The UV magnitude at constant number density M{jy{z) 
brightens with decreasing redshift o ver this redshift rang e 
(similar to the evolution in M{jv, see lBouwens et al.ll2010d ). 
The values of Muv(z) are based on the integral of the lu- 
minosity functions, and they are more robust than the indi- 
vidual luminosity function parameters. While there is non- 
negligible uncertainty on the parameters of the UV lumi- 
nosity functions in the literature, this has little effect on 
the integral of the luminosit y function. Using the z = 8 UV 
luminosity function from Mc Lure et al.l (|2010f ) gives an ab- 
solute magnitude M{jy{z) within 0.1 mag of the value in 
table [1] 

Figure [T] illustrates an importance of using galaxy sam- 
ples on the basis of constant number density rather than 
constant luminosity. Objects at constant luminosity (or stel- 
lar mass) correspond to very different number densities at 
different redshifts. The evolution in the UV luminosity func- 
tions corresponds to a change in number density of about a 
factor of 5 from 2 = 7 to z = 3 at fixed absolute magnitude, 
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Table 1. Properties of galaxies with comoving number density n = 2 X 10 * Mpc 



M 



* b 



ref.' 



^i 



ref/ 



M' 



log If _ 



'M/ 



ref.' 



Original IMF£_ 



Derived from the UV luminosity functions 



2.3 


-20.70 


-1.73 


2.8 


RS09 


1.71 


B09 


-21.06 1.64 


9.9 


R06,Sh05 


Salpeter 


3.1 


-20.97 


-1.73 


1.7 


RS09 


1.77 


B09 


-21.1 1.69 


10.0 


R06,Sh05 


Salpeter 


3.8 


-21.06 


-1.76 


1.1 


B07 


1.67 


B09 


-21.0 1.60 


9.6 


St09 


Salpeter 


5.0 


-20.69 


-1.69 


0.90 


B07 


0.71 


B09 


-20.6 1.03 


9.3 


St09 


Salpeter 


5.9 


-20.29 


-1.77 


1.2 


B07 


0.06 


B09 


-20.3 0.67 


9.0 


St09 


Salpeter 


6.8 


-20.11 


-1.94 


0.90 


BIO 


0.0 


FIO 


-19.9 0.50 


8.7 


St09 


Salpeter 


8.0 


-20.28 


-2.00 


0.38 


BIO 


0.0 


FIO 


-19.4 0.29 


8.6 


FIO 


Salpeter 










Derived from the Stellar 


mass functions 








0.1 


10.96 


-1.18 


3.1 


M09 








11.2 




Kroupa 


1.5 


10.91 


-0.99 


1.0 


M09 








11.0 




Kroupa 


2.5 


10.96 


-1.01 


0.40 


M09 








10.7 




Kroupa 


3.5 


11.38 


-1.39 


0.053 


M09 








10.3 




Kroupa 



'^ Median redshift of each sample. 

Characteristic magnitudes (masses) in units of magnitude (solar masses), faint end slope, and space density at the 
characteristic magnitude (mass) in units of 10~^ Mpc""^ ma g"^ (10~^ Mpc~ ^ dex~ ^). 

'^ References for lumi nosity or mass function para meters: iReddv &: Steidell l|2009l . RS09), iBouwens et al.l l|2007l . B07), 
iBouwens et al.l ll2010al . BlOV lMarchesini et af] ||2009| . M09). 

" Average attenuation at 1600 A at ML^ i n magnitudes. ^_^ ^^^^^_^^^^^^ ^_^ 

^ References for the attenuation at 1600 A: IBouwens et al.l l|2007l . B07'). lFinkelstein et aO l|2010l . FIO). 
J Absolute magnitude of galaxies at comoving number density n{<M!^^) = 2 X 10~* Mpc"'' 
9 Extinction corrected SFR in units of Mq yr~ 



^ derived from M^y 



and Ajjv 



For UV-selected galaxies, these are the mean stellar mass in units of Mq for galaxies with ML^. For stellar-mass-selected 



h 

i„i ... , uv,,^^„v,„ £,u,.u,.„v,„, „.i^uv, u,i^ „,.v, .,,^u,. u„^i.ui ,,!„„„ .11 uim„ „i i„lj i„i j,u,iu,.^iv,„ ..i„ii i"Uy. 1"! U„^ii„i iiiU 

galaxies, these are the stellar masses in units of Mq for galaxies at comoving number density n{>Ml) = 2 X 10"'* Mpc"''. 
^ Reference s for stellar masses for the UV-sel ected galaxies: iReddv et al.l l|2006l . R06), IShaplev et al.l l|2005l . Sh05), 
IStark et ahl | |2009| . St09'). iFinkelstein et al] 1I2OIOI . FIO). All stellar masses have been converted to a Chabrier IMF. 
^ Original IMF for stellar masses. Stellar masses using a lSalpeterl l|l955l ) or lKroupal l|200ll ) IMF were converted to a Chabrier 
IMF by multiplying by a factor of 0.54 and 1.13, respectively. 



Muv = —20 mag. Galaxies selected at constant luminosity 
at different redsliifts will result in samples with very differ- 
ent evolutionary patf is compared to galaxies w ith constant 
number density (e.g.. lConroy fc Wechslerll2009r) . 



2.2 Extinction and star formation rates 

While the rest-frame UV light traces the instantaneous SFR 
of massive stars, it is subject to dust attenuation. Re- 
cent evidence shows that the mean dust attenuation de- 
creases both with decre asing gal axy luminosity and increas- 
ing redshift dReddv et a l. 2008. 2010l : IBouwens et"al]|2009l : 
IFinkelstein et al.ll201Cl ). We quantify the amount of extinc- 
tion for galaxies with M(jv(z) using published values for 
the dust attenuation at rest- frame 1600 A, Auv, for star- 
forming galaxies as a function of luminosity and redshift, 
listed in table [T] As there is little evidence for dust attenua- 
tion at a > 6 l|Bouwens et al.ll2009l . l2010d : IFinkelstein et all 
I2OICI ) we assign ^uv = mag for these galaxies. 

We combine the UV absolute magnitudes and dust at- 
tenuation to derive dust-corrected SFRs, ^{z), 



^lt^\ ^ -|^Q0.4[A[;v-48.6-Afuv{z)-M(z)l 

xA-KdL(zf X 7.8 X 10"^^ MQyr"\ 



(2) 



nosity distance in units of cm. Th e numerical fact or is from 
iKennicutd (|l998i) converted to a lChabrieiJ l|2003l ) IMfQ 

There are several sources of uncertainty on the derived 
SFRs. The scatter in the dust extinction at M(jv is non- 
negligible for galaxies at the lower end of our redshift range 
where the luminosity functions are well measured, producing 
uncertainties of a factor of order 2. There is little dust atten- 
uation at the high end of our redshift range, but the luminos- 
ity function is less- well constrained, leading to uncertainties 
of a factor of order 2. One may expect these uncertainties to 
be lower when averaged over the galaxy population at each 
redshift. Lastly, the uncertainty in the UV luminosity-SFR 
calibration adds an additional 0.1 dex to the error budget 
(see footnote 2). Combining these, we apply a scatter of 0.3 
dex to all SFR values for galaxies at 3 < 2; < 8 in this study. 



where ^{z) is the distance modulus, and dL{z) is the lumi- 



^ Thro ughout we deriv e SFRs and stellar masses assuming the 
IMF of lChabrieJ | |2003|) . The stellar mass and SFR of a stellar 
population with a Chabrier IMF are factors of 1.8 lower than for 
the IM F oflSalpetej l ll955^ . and factors of 1.13 higher than for 
IMF of lKroupalll200lD . We note that calibration between the SFR 
and UV luminosity depends on the assumed S FH. We adopt th e 
relation above, which assumes a constant SFH llKennicutt^ll998a^ . 
The SFR calibration for the UV luminosity for the SFH we derive 
below (equation Is} is within 30% for the stellar population ages 
considered here. 



© 2010 RAS, MNRAS OOO.nHTSl 



Rising Star- Formation Histories and Gas Accretion of Distant Galaxies 5 



2.3 Stellar masses 



Various studies have derived stellar masses, M*, for high- 
redshift ga laxies by modelinR the rest-frame UV-to-near-IR 
data (e.g., 'Shapley et al."2005'; 'Papovich et a ll l2005l. [20061 : 
Rcddv ct al. 2006; Stark ct al. 2009; Finkclstci n"et al.ll2010l : 
Labbe et al.ll2010al ). We use these results to derive the me- 



dian stellar mass for galaxies at each redshift brighter than 
M{jy{z). We present these stellar masses in table [l] Both 
the scatter in the stellar mass distribution at M{jy and 
systematics in the stella r mass estimation contribute to 
the e r ror budget (see e.g.. | Papovich et al.l 120061 : IStark et al.l 
I2OO9I : iGonzalez et al.l bOloD . At z < 6 the scatter in the 
SFR-stellar mass relation dominates the ran dom errors, and 
we a dopt the scatter dlog(M*) — 0.5 dex (JGonzalez et al.l 
[2OI0I). At 2 > 6, random errors from the stellar-population 
modeling are comparable to the scatter in the SFR-stella r 
mass relation (JFinkelstein et al.ll2010l : lGonzalez et al.ll2010f ). 
At z ^ 7 and z ^ 8 we ado pt a scatter of dlog(M*) = 
0.5 dex, (JGonzalez et al.ll20ld l. However, at z = 7 the 68% 
the lower bound on the confidence range extends to a fac- 
tor of order 5 and at z ~ 8 it extends to a factor of order 
10 iJFinkelstein et al.ll20ld l. We adopt these values for the 
scatter here. 



We also derived the stellar mass for galaxies with n(> 
Mi) = 2x 10~* Mpc" -^ at < z < 4 by i ntegr ating the stel- 
lar mass functions of iMarchesini et al.l (|2009l ) in the same 
way as for the luminosity functions above. Table[T]lists these 
stellar masses, Mi{z), derived from the mass functions, and 
converted to a Chabrier IMF. At z < 3 the parameters of 
the stellar mass functions are well constrained for the stellar 
masses of interest here, and the formal uncertainty on the 
stellar mass at fixed number density derived by integrating 
the luminosity function is small, <0.1 dex. At 3 < z < 4 the 
stellar mass function is less well constrained, and we find 
that the uncertainty on the mass function corresponds to 
uncertainties on the stellar mass at n{> M") = 2 x 10~ 
Mpc"^ of log(M^/M0) = 10.2^0:7. In addition, we adopt 
a min imum error bar of 0.3 dex (see, e.g., Marchcs ini et al.l 
I2OO9I ) on the stellar masses at fixed number density calcu- 
lated by integrating the stellar mass functions. 



The stellar masses derived at constant number density 
from the mass functions provide an important comparison 
sample for the UV-selected galaxies. At z <3, observations 
suggest that there exist both infrar ed luminous and UV - 
luminous star-forming galaxies (e.g.. IPapovich et al.ll2006h . 
which make a substantial contribution to the cosm ic SFR 
history (;e.g.. lReddv et al.ll2008l : lBouwens et"alll2009l ). There 
also exist a.t z < 3 galaxies with suppressed or absent 
star-formation (e.g., Kriek et al.l 120061 : Ivan Dokkum et al.l 
l2009l : IPapovich et al.ll2010l ), although such systems are much 
rarer or absent a t z > 3.5 IjBrammer fc van Dokkumll2007l : 
IStutz et al]l2008l ). The UV luminosity functions at z < 3 
may be incomplete for these types of galaxies. The stel- 
lar masses at constant number density derived from the 
IMarchesini et al.l ([20091) mass functions provide a separate 
measure of the evolution of the stellar mass growth in galax- 
ies selected by means other than the UV luminosity. 
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Figure 2. The SFR of galaxies at constant number density n = 
2 X 10"'* Mpc~^ as a function of lookback time measured from z = 
11 to the observed redshift, z. The SFRs are derived from the UV 
luminosity and corrected for dust extinction as discussed in 5 12.21 
At this number density the SFR increases with time (decreasing 
redshift) from 2 = 8 to z = 3, with data points indicated as filled 
symbols. The red solid lines and shaded region show fits 'l'(i) ~ t" 
to the data from 3 < z < 8 with parameters and 68% confidence 
range given in the text. The power-law fit formally has a better 
goodness-of-fit than the black dashed line, which shows a model 
where the SFR increases exponentially in time, ^{t) ~ exp(4/r). 
The SFR unambiguously increases smoothly with time from 2 = 8 
to z = 3. Galaxies at z < 3 appear to depart from this evolution, 
and we indicate the datum at z = 2.3 with an open symbol. 



3 THE STAR FORMATION HISTORIES OF 
DISTANT GALAXIES 

Figure [2] shows the evolution of the dust-corrected SFR av- 
eraged over galaxy populations with n — 2x 10~* as a func- 
tion of lookback time, calculated from a formation redshift 
Zf — 11 to the observed redshift jj using the data in table [l] 
There is an unambiguous smooth increase in the SFR from 
z = 8 to 3. The relationship is consistent with a power law, 
described by 



*W 



IM, 



yr" 



{t/rr 



(3) 



with a = 1.7 ± 0.2 and r = 180 ± 40 Myr fitted over the 
data from z = 8 to 3. We adopt here the power-law fit as 
it formally has a better goodness-of-fit compared to other 
simple functions, although it is not a unique representation 
of the data. Fitting a model to the data where the SFR 
increases exponentially in time, 'I'(z) = 5'oexp(f/r), gives a 
minimum x^ ~ 1-1 for 4 degrees of freedom (for r = 420 Myr 
and ^0 ~ 1-4 Mq yr~^), compared to x^ ~ 0.66 for 4 degrees 
of freedom for the power-law fit above. However, statistically 
these fits are indistinguishable statistically. 



•^ The choice of formation redshift Zf = 11 is arbitrary and does 
not afi'ect our conclusions. Changing the formation redshift adds 
only a constant offset in time to equation [3] 
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Here, we concentrate on the evolution of galaxy prop- 
erties at constant number density n — 2x 10"* because this 
describes the evolution of "typical" galaxies (§ 2), and this 
is where the UV data are most complete (see references in 
table [l]). However, we also studied how the SFH varies with 
galaxy number density, ranging from n — 1 x 10~^ Mpc~^ 
(corresponding to fainter galaxies with ~Z/* + 1 mag) to 
n = 4 X 10~^ Mpc~^ (corresponding to brighter galax- 
ies with ~L* — 0.5 mag). We observe a modest change in 
the slope of equation [3] for these samples, from a = 1.5 at 
n = 1 X 10"^ Mpc"^ to a = 1.9 at n = 4 X 10"^ Mpc"^, 
but these are consistent within the measured uncertainty, 
5a=Q.2. However, the normalization (t~") changes relative 
to that for n = 2 X 10"" Mpc"^ by a factor of 1/22 at 
n = 1x10"^ Mpc"^ to a factor of 12 at n = 4x10"^ Mpc~^. 
Therefore, over the range of number densities considered 
here (n = 4 x 10~^ — 1 x 10~^ Mpc"^) galaxies experience 
SFRs that increase as a power-law in time with a slope, a, 
that is approximately constant. The SFR scales only by a 
multiplicative factor related t o number density (o r mass, 
similar to the simulations of iFinlator et al.l |2010|) , where 
galaxies with higher number den sity (lower UV lumino sity) 
experience slower evolution (e.g. jBouwens et al.ll2010al ). 

At z < 3 the SFRs derived for UV-selected galaxies ap- 
pear to depart from the SFH in equation [3] (figure [2]). At 
2 ~ 2 — 3 galaxies with UV luminosities near L* show a sig- 
nificant heterogeneity in their dust properties, and galaxies 
at fixed UV luminosity include both intrinsically less lu- 
minous galaxies with less dust extinction and intrinsically 
more luminous gal axies whose UV emission is obscured by 
dust IjReddv et al.l.2010. ). and so we may expect the depar- 
ture from the derived SFH. However, our derived SFH is 
consistent with the space density of Ji'-selected ultralumi- 
nous IR galaxies (ULIRGs) at 1.5 < z < 2.5, which have 
n = 2 x 10"" Mpc"^ and (*) = 120 Mq yr"^ (|Daddi et all 
120071 ). Therefore, while UV-selected samples seem incom- 
plete for star-forming galaxies at redshifts z < ?>, the de- 
rived SFH remains consistent when we include observations 
of obscured galaxies. 



Therefore, galaxies with this constant number density 
have average SFRs that increase with time for redshifts 
3 < z < 8. This is sim i lar to expectations from theory (e.g., 
Robertson et al.l |2004|; iDe Lucia et all l2006l: IFinlator etahl 



0.5 1 



Redshift z 
2 3 4 5 6 



20061 . I2OIOI : iBrooks et ai]|2009D . which predict rising SFHs 



over this redshift range. We stress that the SFHs we derived 
here describe the cosmologically average evolution of the 
galaxy population at these redshifts, and individual galaxies 
likely experience unique SFHs, including stochastic bursts 
and periods of reduced star-formation. In particular, at later 
times, z < 3, the SFHs of galaxies seemingly plateau, and 
at later times galaxies must experience SFHs that decline 
with time in order to explain the existence of lower redshift 
galaxies with evolved stellar populations evidenced by the 
existence of stron g 4000 A/Balmer breaks by 2: ~ 2 (e.g., 
iKriek et al.l 120061 ). Nevertheless, as we discuss below (§|31), 
the average smoothly rising SFHs for the galaxy population 
at early times reproduce the stellar mass growth in galax- 
ies. We therefore, conclude that galaxies with this constant 
number density experience rising SFHs during the epoch de- 
marcated by 3 < 2 < 8. 
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Figure 3. The stellar mass of galaxies as a function of lookback 
time for galaxies with constant number density. Filled squares 
show the stellar masses for galaxies with UV luminosities such 
that n(< A/(jv) = 2 X 10~" Mpc-^ from x = 8 to 3. The 
heavy red line and shaded region shows the possible range of 
stellar mass evolution expected from SFH measured in figure (2] 
This line is not a fit, as we have made no adjustments to the 
normalization or slope of this line. The stellar mass growth ex- 
pected from the measured SFH agrees with the data for UV- 
selected galaxies at 3 < 2 < 8. The measured masses for UV 
samples at 2 < 4 are slightly below the expected evolution for 
possible reasons discussed in the text. The filled circles show the 
evolving stellar masses for g alaxies with stellar mas s such that 
n{>M!,) = 2 x 10"" Mpc-3 JMarchesini et al.ll2009t ). These fol- 
low the expected stellar masses from the SFH for redshifts 2 > 1. 
Therefore, the SFH and stellar mass growth are consistent over 
this redshift range. 



4 THE EVOLVING SFR-STELLAR MASS 
RELATION 

The measured SFH provides a direct prediction for the stel- 
lar mass evolution for galaxies at this constant number den- 
sity. For a SFR that evolves as \&(i) ~ i" (equation (3)) the 
stellar mass evolves roughly as Mt{t) ~ J 'if(t)dt ~ t'"+^'. 
However, one must account for stellar-mass loss due to bary- 
onic gas returned to the ISM of a galaxy from stellar evolu- 
tion. We therefore computed the expected stellar mass evo- 
lution using_tlw_SFH_^wn equation [3] with the 2007 version 
of the lBruzual fc CharlotI l|2003r ) stellar-synthesis model. We 
use an input stellar population with Z = 0.2 Zq metallicity 
and a Chabrier IMF with no chemical evolution. 

Figure |3] shows the evolving stellar mass, M^, as a func- 
tion of lookback time for galaxies with n = 2 x 10"" Mpc"''. 
The measured masses of the UV-selected galaxies at this 
number density are shown as black squares. The heavy red 
line and shaded region in figure |3] shows the expected increase 
in the stellar mass from the range of SFHs in figure [51 It is 
important to note that the red curves in figure [3] are not fits 
to the data, but rather predictions from the derived SFH. 
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The stellar mass evolution of the UV-selected galaxies 
at 3 < 2 < 8 agree well with this expected growth from 
the measured SFR history. At z ~ 8 the stellar mass for 
UV-selected galaxies with this number density lie above the 
expected value, although still within the scatter. 

At 2 < 4 there is slight evidence that stellar masses 
for the UV samples fall below the expected trend by about 
0.3 dex (although again within the scatter). This could be 
symptomatic of older stellar populations with high mass-to- 
light ratios which are plausible to the limits o f the data 
(e.g., iPapovich et"ai] 1200 ll : IShaplev et all |2005| ). Alterna- 
tively, the stellar masses derived for the UV samples as- 
sumed SFHs that are either constant or decline with time, 
rather than increase with time. As discussed in Appendix iBl 
stellar masses from rising SFHs are larger by «0.2 dex for 
galaxies at z ~ 3 — 4 compared to mode ls with SFRs tha t 
are constant or decline with time (see also iLee et al.ll2010br ) . 
In contrast the stellar masses from models with rising SFRs 
are only <0.1 dex higher at z ~ 5 — 6. These shifts are indi- 
cated by the open boxes in figure O and they would account 
for most of the observed offset. Furthermore, this offset may 
be also indicative of the fact that at z ~ 2 — 3 galaxies with 
UV luminosities near L* show a sign ificant heterogeneit y in 
their dust properties (see § O and iReddv et al.ll201(j ). or 
as a result of the scatter in the UV-luminosity-dark mat- 
ter relation (see Appendix [X|) . Lastly, by z < 3 there exist 
mass ive galaxies with "suppressed" SFRs (e.g.. lKriek et al.l 
120061 ). which could be missed by the UV selection. These 
effects may bias the derived stellar masses for these galaxies 
to lower values. In contrast, these effects are less significant 
at higher redshifts z > 4 because there has been less time 
for them to manifest. 

Because of potential biases in the UV-luminosity- 
selected sample, figure [3] also shows the stellar masses for 
galaxies at z < 4 with constant number density, n(> 
Mj) = 2 X 10~^ Mpc~'^ , derived from the mass functions of 
iMarchesini et al.l (|2009l ). The selection based on stellar mass 
mitigates most of the biases for the UV galaxies. Remark- 
ably, the evolution of stellar mass for galaxies at constant 
number density agrees well with the stellar mass growth ex- 
pected from the SFH for z > 2. This supports our conclusion 
that the cosmologically averaged SFHs of distant galaxies 
are smoothing rising from z = 8 to 3 because this matches 
both the increased UV luminosities and stellar masses of 
galaxies with this number density. 

The models with rising SFHs contrast assump- 
tions that distant galaxies experience SFRs that 
are either constant or d e creasi n g with time (e.g. , 
Papovich etai] l200ll. [2OO5I. I2OO6I: Ishaplev et al.l 1 2001 : 



rr 



Forster-Schrei ber et al .l 120041: IShaple v et al.l 120051: 
Reddv et alBoOG: Forster-Schreiber et al.r2009: Sta rk et all 



20091 : iFinkelstein et al.ll2010l : iLabbe et al.|[20I0ai ibi). Several 



other studies have suggested SFHs that increase with time 
based on the evolution i n the SFR-stellar mass relation 
at high redshifts (z < 6) l|Renzinil [2OO9I : IStark et al.ll2009l : 
iMaraston et al.l 20ld) and also based on theoretical motiva- 
tion ( Finlator et al.ll2007l .l2010l. Our analysis demonstrates 
unambiguously that the cosmologically averaged SFRs for 
typical galaxies increase smoothly with time (decreasing 
redshift). Individual galaxies likely deviate from this 
"smoothly" increasing SFRs because they have forma- 
tion histories that involve stochastic events with discrete 



changes in the instantaneous SFR (e.g., Baugh et al.l boosi : 
iFinlator et aLll2006l . l2007l : ISomerville et al.ll2008l ')."and they 
may experience SFHs that deviate from the average history 
derived here. The SFHs that we derive here corresponds to 
a cosmologically averaged evolution that accounts for the 
observed evolution of the galaxy populations. 

The self-consistency between the rising SFHs and the 
stellar mass growth constrains jointly (degenerately) both 
the star formation duty cycle and the slope of the high- 
mass end of the IMF. Observations of clustering measure- 
ments and the correlation and scatter of the SFR-stellar 
mass relation of high-reds hift galaxies imply a high star- 
formation duty cycle ( e.g., iDaddi et al.l 120071 : iLabbe et al.l 
l2010al : iLee et al.ll2010al ). If the duty cycle is near unity for 
galaxies at 3 < z < 8, then the consistency between the 
rising SFHs and stellar mass growth here is evidence that 
the high-mass end of the IM F in these galaxie s mu st have a 
slope consistent with that of IChabrien (12003) and ISalpeten 
(I9551) when averaged over galaxy populations. Otherwise, 
there would be a larger offset in the SFR-stellar mass rela- 
tion, which is not evident in the data. 

If the star formation duty cycle is lower as suggested by 
clustering analy ses of (predominantly) suh-L* sources (e.g., 
iLee et abl |2009| . although see Finlator et al. 2010) and b y 
some models of galaxy formation (e.g.. iBaugh et aLlbOOa ). 
then our integrated SFH requires an IMF weighted toward 
higher mass stars in order to match the measured stellar 
mass growth. For example, a duty cycle of 50% would re- 
quire a change of ~ 0.2 to the slope of the IMF to produce a 
higher ratio of the number of 0-type stars to the number of 
G-type stars compared to that produ ced by a Salpeter-like 
IMF. This is similar to the findings of lBaldrv fc Glazebrookl 
(l2003t ) based on SFHs derived from the UV to IR luminosity 
densities of low-redshift galaxies. Such top-heavy IMFs have 
been also proposed to explain the joint cosmic evolution of 
the stellar mass and SF R densities fe.g.. lBaugh et al.ll2005l : 
lDavg|2008l : iLacev et al.ll200S ). although these are more ex- 
treme IMFs than the change described above. Lastly, our 
results are unable to constrain the low-mass end of the IMF 
as any change in the shape of lower mass limit simply scale 
both the SFR and stellar mass in approximately the same 
way, shifting both the datapoints and the curve in figure [3] 
by approximately the same amount. 



5 IMPLICATIONS FOR THE GAS MASS AND 
GAS ACCRETION RATE 

The evolution in figures [2] and [3] leads to the conclusion that 
galaxies increase both their SFRs and stellar masses over 
time. This simply is not possible if galaxies form with a fixed 
initial cold baryonic gas mass. In this case the SFR would 
decline with time as the gas mass is converted to stars. 

In this section we use our measured evolution in the 
SFRs and stellar masses of galaxies to predict the evolution 
of the gas mass in galaxies with with n = 2 x 10"* Mpc~^ 
at high redshifts. Throughout this section we assume that 
the local gas-su rface density (Sgaa) - SFR surfa ce density 
(Esfr) relation (jSchmidtHigSgl : lKennicuttlll998bl ') applies to 
high-redshift galaxies. 



SsFR 



Mq yr~ikpc 



= A 



Mq pC-2 



(4) 
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where A = 1.7 x 10 * and A'' — 1.4 have their normal 
values (t hough we use the ca hbration for a Chabrier-like 
IMF, see lSomervihe et aLllJOOJ . see also Button et al. 2010). 
The assumption that star-formation in high-redshift galax- 
ies follows this relation remains untested directly at these 
redshifts, and such tests will require direct measurements 
of the HI gas, which is currently in f easible . However, recent 
studies for example bv lDaddi et al.l (|2010a ) suggest that the 
star-forming efficiency in normal star-forming disk galaxies 
is consistent with this assumption, though it may not apply 
in extreme classes of star-forming galaxies. As the "typical" 
galaxies studied have SFRs and properties consistent with 
"normal" star-forming disk galaxies at 2: ~ 2, our assump- 
tion is reasonable. 

We multiply equation [4] by the areal size of the star- 
forming region in galaxies, approximated by a = Trr^/j; 
where ri/2 is the familiar half-light radius, which relates the 
gas-mass to the SFR, 



*(2) 



1 M© yr- 



C 



ri/2{z) 



Ikpc 



A'/gaa(2) 



lO^Mp 



(5) 



where the constant of proportionality is C = 1.7. While the 
statistical uncertainty on equation [5] is 0.12 dex (jKennicutg 
[l998bj), the true uncertainty may be larger at high redshifts 
where these relations have not been measured directly. 

The rest-frame UV size ev olution of high-redshift galax- 
ies has been characterised bv iFerguson et al.l (|200J), who 
showed that the size evolution for star-forming galaxies from 
2 ~ 5 to z ~ 2 is consistent with the growth of galaxy disks 
that scale with the size of their dark matter haloes ( Mo ct ah! 
119981 ). r{z) oc H{z)~^, where H(z) is the Hubbl e constant 
at redshift z (consistent with lOesch et al.ll2010al ). We nor- 
malize this relation to produce the characteristic half-light 
size at z — 4: from Ferguson et al., r(z=4) — 1.7 kpc with 
a scatter of approximately 0.6 kpc. The galaxy sizes then 
evolve as 



r{z) 



1.7 



H{z) 



H{z = 'i) 



kpc. 



(6) 



where H{z—4) = 430 km s ^ Mpc ^ . The extrapolated size 
evolution from this equation reproduces the sizes measured 
for z = 7 galaxies in HST /WFC3 data, r(z=7) = 0.85 kpc, 
consistent with the observations of lOesch et al.l (|2010al ). 

Combining equations [5] and [S] and solving for the gas 
mass yields 



MgaB(z) 

5.9 X 10« Mq 



*W 



i5/7 r 



IM, 



©yr" 



r(z=4) 



-i4/7 



1.7kpc 

H{z) 



(7) 



H{z=4) 



-4/7 



The uncertainty on the gas mass includes the scatter in the 
galaxy size measurements, and includes the uncertainty from 
the SFR-gas mass surface density relationship (a factor of 
(JlogAfgas ~ 0.11 dex). However, as discussed above, equa- 
tion [7] implicitly assumes the star-formation law from the 
relation given in equation[4]and neglects any uncertainty on 
this relation at high redshift. 
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Figure 4. Predicted redshift evolution of the gas accretion rate, 
Afacc and the gas mass, the derived stellar mass inferred from the 
empirical SFH, and the empirically derived evolution of the SFR 
for galaxies with constant number density n = 2 X 10~ Mpc~ . 
Here, Mace is the net gas accretion rate and does not distin- 
guish between the contributions of gas inflow and outflow. The 
gas mass and gas accretion rate are discussed in the text, and 
assume the Kennicutt— Schmidt conversion between gas and star- 
formation surface density. Each curve shows an evolving quantity 
as indicated in the legend. Gray-shaded regions show redshifts 
where there is no data, and the curves are extrapolations. The 
left abscissa shows the scale for the gas accretion rate, the SFR, 
and the rate of gas mass returned from stellar remnants. The right 
abscissa shows the scale for gas mass and stellar mass. The down- 
ward arrow indicates the redshift where the stellar mass equals 
the gas mass. The upward arrow indicates the redshift where the 
SFR equals the gas-mass accretion rate. 



5.1 Gas mass and gas accretion rate 

Figure |4] shows the derived redshift evolution of the stellar 
mass (Mi) inferred from the empirical SFH (§|4]), and the 
predicted redshift evolution for the gas mass (Mgag) from 
equation [7] for galaxies with constant number density n — 
2 X 10"'* Mpc"''. The stellar mass increases with time, as 
discussed in § |3] The gas mass also increases from z = 8 
to z = 3, which is required to fuel the increasing SFRs. 
Our empirical relations predict that from z = 10 to z = 
8 the gas mass is an order of magnitude larger than the 
stellar mass for galaxies with n = 2 x 10~* Mpc"''. The 
gas mass continues to increase at z < 8. While the stellar 
mass initially lags behind the gas mass, it increases quickly 
in response to the SFR. Around 2 »; 4, the stellar mass 
equals the gas mass for galaxies with this constant number 
density, and extrapolating to lower redshift the stellar mass 
becomes the dominant baryonic component for galaxies at 
this number density. 

We use the derivative of the gas mass and the SFR 

ace, for these galax- 

be driven by either 

smooth accretion (e.g.. iKeres et al.l 20091), or gas delivered 



to predict the gas accretion rate, M; 
ies. The gas accretion rate could 



via significant mer gers (e.g., [Robertson fc Bullockl I2OO8I : 



IStewart et al.l 120081 ). and differences in these processes do 
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not affect the simple empirical relations we derive here. Here 
we do not difTerentiate between possible physical processes 
that contribute to the net gas accretion rate but rather we 
define this to be simply the net rate at which galaxies ac- 
quire their baryonic gas. Therefore, the net gas accretion 
rate here includ es the effects of both gas infall and ga s out- 
flows (see, e.g.. iDave erai]|2006l : [button et al.lboiol '). The 
time-rate-of-change of the gas mass has several components. 



Mgas(i) = Macc(i) - *(f) + Afrcc(t), 



(8) 



where M^ec is the amount of gas returned ( "recycled" ) to the 
galaxy from stellar remnants (taken from the 2007 version 
of the Bruzual & Chariot 2003 models). 

Figure U] shows the predicted gas accretion rate for 
galaxies with number density, n — 2 x 10~^ Mpc""^ using 
equation[8]. The figure also shows the SFR and gas-recycling 
rates. The gas accretion rate increases from z = 8 to 3 and 
keeps pace with or exceeds the measured SFR. At these red- 
shifts these galaxies acquire gas at least as rapidly as it can 
be converted into stars - this is the "gas accretion epoch" 
for these galaxies, and is similar to inferences drawn from 
SPH simulations of galaxy formation (Keres et al. 2005). At 
2 < 4 the gas-accretion rate declines, and galaxies consume 
their gas at a rate faster than they acquire it. 

Extrapolating the trends in figure |4] to lower redshifts 
{z < 3), our model shows that the stellar mass of gala^c- 
ies with n — 2 x 10~* Mpc~'^ exceeds the available gas 
mass and the SFR exceeds the gas-accretion rate. These 
trends can not continue ad infinitum as the available gas 
will eventually be consumed. Indeed, figure [3] indicates that 
at redshifts z < 2 — 3, the stellar masses of galaxies with 
n — 2 X 10~^ Mpc~^ depart from the expected growth 
based on the evolving SFR, and this appears consistent 
wit h the evolution of neutral gas in absorption- line stud- 
ies (jProchaska fc Wolfe 2009). By 2: ~ 2 some galaxies with 
the nu mber density cons idered here have low or suppressed 
SFRs (jKriek et al.ll2006l ). which presumably have low cold- 
gas fractions. This implies that the gas accretion rate is de- 
creasing at these redshifts, consistent with galaxy forma- 
tion simulations (Keres et al. 2005,). Several physical effects 
and feedback processes have been propo sed to prohibit gas 
cooling and accretion in gal a xies (e.g., Birnboim fc Dekej 

' ' 12'ooi 



20031: ISpringel et al.l l2005al: iHopkins et al 



2008 



Croton et al.ll2006l : ISomerville et al.ll2008l l. which should ap 



ply at the mass scales of galaxies with this number density 
at these redshifts. 

The semi-anal ytic galaxy formation models of 
iDutton et al.l (|2010|) find that star-forming galaxies at 
high redshift have SFRs in a steady state with the total gas 
mass rate-of-change (gas accretion minus gas losses due to 
outflows) of Afgas ~ ^(i)- The gas accretion rate predicted 
by our model seems consistent with this conclusion (al- 
though large disparities between the g a s accr etion rate and 
SFR are not excluded). iDutton et al.l (|2010|) interpret the 
relative consistency between the SFR and net gas accretion 
rate as a result of increased gas densities in galaxies. Our 
empirical results may support this scenario. 



5.2 Predictions for galaxy gas-mass fractions 

We determine the gas fraction as /gas ~ A/gas (2) x [Afgas(2)-|- 
Mt.{z)]~^ using the evolution of the stellar mass inferred 
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Figure 5. Empirical constraints on the gas fraction, /gas = 
Afgas X (Mgas + Mt)~^. The hatched region shows the expected 
range of average gas fractions for galaxies with constant num- 
ber density n = 2 X 10~* for z = 8 to ^ = 3. The dashed 
lines indicate these relations extrapolated to z = 1.5. Data 
points are labell e d, and include UV-select ed galaxies (filled boxes, 
lErb et al.l I2OO6I: iFinkelstein et al.l [20 09') and sub-mm-selected 
galaxies (filled ci rcles, iKneib et al. 2005; Papovicl^_e^^lJ_[2009|; 
IPaddi et a .Il2009[). The data points for iForster-Schreiber et al.l 
1120091 ) and Tacconi et al.l yOlO) show the mean value for galaxies 
with 2^2 and M« > 3 X lO'^ Mq, the approximate stellar mass 
atz ^ 2-2.5 such that n(>M») = 2xl0~* Mpc'^. Our empirical 
relation generally accounts for the measured gas fractions of these 
galaxies. There is evidence of mass-dependent evolution. Galax- 
ies with higher stellar masses have lower gas fractions (e.g., the 8 
o'clock arc, Finkelstein et al. 2009; GN20, Daddi et al. 2009), im- 
plying faster gas-consumption. Likewise, lower mass galaxies have 
higher gas fractions than what is predicted from the extrapolated 
relationship for galaxies with n = 2 X 10"** Mpc~^ 1 Daddi et al.l 
2010a, see also the lower mass samples of Forster-Schreiber et al. 
2009 and Tacconi et al. 2010) . 



from the SFH (§ |4)| and the evolution of the gas mass from 
equation [7] The hatched region in figure [S] shows the pre- 
dicted range of the gas fraction from 2 = 8 to 2 = 3 for 
galaxies with constant number density n — 2x 10^^ Mpc"''. 
The dashed lines show the relation extrapolated to lower 
redshift. The gas mass fraction for these galaxies decreases 
with decreasing redshift, parametrised as /gas oc (1 -I- z)^ 
with F ~ 0.9. This follows from the evolving stellar mass 
and gas mass derived above. At 2 ~ 3 — 5 the range of possi- 
ble gas-mass fractions reach 50%, similar to our conclusions 
from figure m 

Figure [S] shows that the predicted gas fractions 
are consistent with available observations of galaxies at 
2 '^ 2 with stellar masses such that n(> Mi) = 2 x 



10" 
from 



Mpc , including galaxies with gas masses derived 



C O measurement s Kn^bet_aJj_j2005|); [Dad^i_et_aJj 



20091) : iPapovich et"all (|2009l ): iTacconi et al.l (|201G| ) and 



those derived from the SFR surface density inferred from 
Ha measurements (lErb et al.l I2OO6I : IFinkelstein et al.l l2009l : 
iForster-Schreiber et al.l [2009|) . There is evidence that the 
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evolution of the gas fractions is mass dependent. Figure [5] 
also shows gas fractions of several galaxies with masses much 
greater than the galaxies with n — 2x 10~ Mpc~ (these are 
presumably rarer systems, with lower number density) . This 
includes the gravitationally lensed UV-selected, "8 o'clock 
arc" at 2 = 2.7 with high stellar mass (> lO'^^ Mq) and low 
gas fra ction, which may have experienced accelerated evo- 
lution (Finkolstein et al. 2009'). Similarly, the massive sub- 
mm- sclccted galaxy GN 20 at z — 4 has a low gas-mass frac - 
tion JDaddi et al.1 120091 ') (see also. iMichalowski et ahlboiol ). 
These systems show evidence of mergers, enhancing the 
;as-consumption ra te and/or the star-formation efficiency 
Daddi et al.l l2010b). Because the empirical gas fractions we 
derive for galaxies with n = 2 x 10~* Mpc"'' agree reason- 
ably with the available observations, galaxies experiencing 
enhanced gas-consumption rates because of significant merg- 
ers seem to be the exception, not the norm at these redshifts. 
Galaxies with lower stellar masses than those selected 
with constant number density, n = 2 x 10~'^ Mpc~^, ap- 
pear to have higher gas-mass fractions l|Daddi et al.ll2010al . 
see also the lower mass galaxies in the samples of Forster- 
Schreiber et al. 2009 and Tacconi et al. 2010). This is further 
evidence suggesting that the gas-accretion and gas consump- 
tion time-scales dependent on galaxy mass. 

Our empirical results make very clear predictions for the 
cosmologically averaged gas masses of galaxies with number 
density n = 2 x lO"'' Mpc"'', and also that the gas fractions 
should increase as function of redshift. While the predicted 
gas masses are consistent with the limited available data, 
more observations are needed. Future observations from the 
Atacama Large MilUmeter Array (ALMA) will have the sen- 
sitivity to detect emission from the predicted gas reservoirs 
in these galaxies, which will test the predictions on the gas 
masses we make here. 



stochastic events and unique SFHs that deviate from the 
average evolution derived here. 

We show that the measured stellar masses for galax- 
ies selected with this constant number density are consis- 
tent with the expected stellar mass growth from the derived 
SFH. If star-forming galaxies at high redshifts have star- 
formation duty cycles near unity, then this implies that the 
high-mass end of the stellar initial mass function is approxi- 
mately Salpeter. For lower duty cycles, our results would re- 
quire an IMF with high-mass slope weighted toward higher 
mass stars (by « 0.2 for a duty cycle of «50%). Otherwise, 
there would be a larger offset in the SFR-stellar mass rela- 
tion, which is not evident in the data. 

We interpret the relation between SFR and stellar mass 
as a result of net gas accretion (gas outflows combined with 
gas infall from either smooth accretion or delivered by merg- 
ers) coupled with the growth of galaxy disks, where the SFR 
depends on the local gas-surface density relation. This leads 
to predictions for the galaxy gas masses and gas accretion 
rates on galaxies z = 8 to z = 3. The gas fraction evolves at 
a rate of /gas oc {l + zf , with F ~ 0.9, and we show these gas 
fractions are in reasonable agreement with available data at 
z ~ 2. These relations imply that at z > 4 galaxies selected 
at this constant number density acquire baryonic gas at least 
as fast as it can be converted into stars. This is the "gas ac- 
cretion epoch" for these galaxies. Interpreting this in the 
context of a semi-analytic model (jDutton et al.l 2010) this 
suggests that the smoothly rising SFHs are a consequence 
of increased gas densities in these galaxies. At z < 4 the SFR 
in these galaxies overtakes the gas-accretion rate, indicating 
a period where these galaxies consume gas faster than it is 
acquired. At redshifts z < 3, galaxies depart from these re- 
lations implying that gas accretion is slowed or prevented at 
later times. 



6 SUMMARY 

We use observed relations of the SFRs and stellar masses 
for galaxies to constrain empirically the SFHs of galaxies 
at 3 < z < 8. We compare the evolution of galaxies at 
these redshifts at constant comoving number density, n(< 
Muv) = 2 X 10"'* Mpc"^. We show that this allows us to 
track the evolution of stellar mass and star formation in the 
direct predecessors and descendants of these galaxies in a 
relatively meaningful way that is not possible using samples 
selected either at constant luminosity or constant mass. We 
summarize our findings as the following. 

The galaxies experience cosmologically averaged SFHs 
with SFRs that grow with time (decreasing redshift) at a 
rate best described by *(t) ~ t" with a = 1.7 ± 0.2. This 
excludes model SFHs where the SFR is either constant or de- 
clines exponentially in time. We find evidence that galaxies 
over a range of number density from n = 4 x 10~^ Mpc~'^ to 
1 x 10~^ Mpc"'^ (corresponding to galaxies with luminosities 
< L\jY — 0.5 mag and < -Luv + 1 niag, respectively), experi- 
ence SFHs with similar power-law increases with time, but 
differ by a multiplicative scale factor. These results are sim - 
ilar to expectations from theory (e.g.. iFinlator et al.ll2010h . 
We reemphasize that these SFHs correspond to the cos- 
mologically averaged evolution of the galaxy populations 
with these number densities. Individual galaxies experience 
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APPENDIX A: ON GALAXIES AT CONSTANT 
NUMBER DENSITY 

In this paper we have studied the evolution of galaxies from 
3 < 2 < 8 at fixed number density. This has the advantage 
that in principle it traces the evolution of the progenitors 
and descendants of galaxies as a function of redshift in a 
way that is not possible using other methods such as selec- 
tion at constant luminosity or mass, quantities which them- 
selves evolve with time. It is also relatively straightforward 
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Figure Al. The evolution of halo ma ss at constant numbe r 
density from the Millennium simulation (jSpringel et al.ll2005a l. 
The datapoints show the halo mass, Mvir, corresponding to con- 
stant number density, n{> M^ir) = 2 X 10~* Mpc~^, at each 
redshift. Halos selected at constant number density grow in mass 
by roughly a factor of 10, similar to the growth in stellar mass 
for galaxies we observe at these redshifts (figure [Sjl- 



to compare results from observational galaxy surve ys to ex- 
pectations from theory. Ivan Dokkum et al.l (|20ld ) applied 
a similar selection to study the size evolution of massive 
galaxies at z < 2. They quantified the effects of mergers on 
their selection using Monte Carlo simulations of the mass 
function, and concluded that selection at constant number 
density produces homogeneous samples at different redshifts 
with low contamination unless haloes grow only via one-to- 
one mergers. Based on theoretical arguements they argued 
that evolution driven by major mergers is highly unlikely 



that evolution driven by major mergers is 
fe.g.. lNaaberai]|2007l : IGuo fc Whit"3l2008l ) 

We studied how well selecting galaxies at constant num- 
ber density tracks the descendants and progenitors of galax- 
ies at different redshifts by exploring the halo-merger trees 
and mass functions from the Millennium simulation (e.g., 
ISpringel et al.ll2005a ). This allows us to study halo growth 
via mergers and accretion in a physically motivated sim- 
ulation that reproduces observational clustering statistics. 
We selected dark- matter haloes with mass, M^ir, such that 
n(>Mvir,2) = 2 X 10"* Mpc"^, for redshifts 3 < z < 7.3. 
Figure lAll shows that the mass corresponding to this con- 
stant number density increases by more than an order of 
magnitude over this redshift range (similar to the growth in 
stellar mass in galaxies we observe in figure |3}- 

We tracked the descendants of the halos selected in 
the Millennium simulation that have halo mass such that 
n(>Mvir,2) = 2 X 10"** Mpc"^ at z < 7.3 and have a 
progenitor at z = 7.3 with mass Mvir such that it has the 
same number density. The top panel of figure IA2I shows the 
"completeness" fraction, which is defined as the ratio of the 
number of the descendants of the galaxies that reside in 
haloes with mass that satisfies n(>Afvir) = 2 x 10"* Mpc"^ 
at lower redshift z compared to the number at 2 = 7.3. 
The completeness fraction is related to the "contamina- 
tion" fraction, which is the fraction of those haloes at red- 
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Figure A2. The top panel shows the completeness of the descen- 
dants and progenitors of dark matter halos at different redshifts. 
The filled squares show the descendant completeness, defined as 
the fraction of haloes selected at redshift z with mass such that 
n(>Mvir) = 2 X 10~* Mpc~^ and a progenitor at z = 7.3 se- 
lected at this number density. The open squares and dashed curve 
shows the completeness for the progenitors of haloes selected at 
constant number density at 2 = 3. The bottom panel shows how 
the completeness changes for galaxies selected at fixed SFR with 
n(>*) = 2 X 10~* Mpc~^. The filled squares show the complete- 
ness for the halos as in the top panel. The dashed cu rve shows 
the completeness for the SAM of JDe Lucia et al.l 1120061 ) using the 
Millennium simulation. The solid curve shows the completeness 
assuming that the scatter in halo mass is d (logMvir) = 0.1 dex 
at fi xed SFR, as pr edicted by the models of JDutton et al.l 1120101 ) 
and JFinlator et~al] jJOig) . 



shift z with mass such that n(>Afvir) = 2 x 10"'' Mpc~^ 
but lacking a progenitor with mass such that n(>Mvir) = 
2 X 10~* Mpc~^ at z — 7.3. Because the number density 
is constant, the completeness and contamination are related 
by (1 — completeness) — contamination. 

The completeness fraction drops from unity at 2 ~ 7.3 
to approximately 0.6-0.7 for 3 < 2: < 6. This implies that 
selecting halos at this number density identifies the major- 
ity of the lower-redshift descendants of objects selected at 
higher redshift. The haloes whose descendants are "miss- 



ing" have masses < 0.5 dex below the mass limit for this 
number density. Arguably if we allowed an evolving num- 
ber density criterion, then we would select more of the halo 
descendants, increasing the completeness fraction. However, 
this would come at the cost of increasing the contamination 
fraction. Similarly, the haloes that contaminate the sam- 
ple have masses at 2 = 7.3 within < 0.5 dex of the mass 
limit for the number density selection. Because these con- 
taminants have masses close to the selection threshold, we 
expect they will have similar evolutionary paths and SFHs 
as the galaxies se lected at constant nu mber density (see § 3 
and discussion in lFinlator et al.ll2010h . 

Figure IA2I also shows the completeness of the pro- 
genitors of haloes selected at 2 = 3 with n(> Mvir) = 
2 X 10"* Mpc~^. The conclusions are similar as those for 
the descendants of the 2 — 7.3 haloes. The progenitor com- 
pleteness becomes lower at higher redshifts compared to the 
descendant completeness, as the progenitors of 2 = 3 haloes 
become slightly less massive than the required value for the 
number density selection. We conclude that selecting halos 
at constant number density traces the majority of the de- 
scendants and progenitors of galaxy halos at 3 < 2 < 7.3, 
and therefore provide samples to study the direct evolution 
in these populations. 

However, thus far we have discussed the completeness 
in the descendants and progenitors of only the dark matter 
haloes. To make the connection to observed galaxies, we 
make the ansatz that there is zero scatter between galaxy 
mass, galaxy luminosity, and halo mass, 



ngal(>I/Uv) = ngal(>Af*) = nhalo(>Mvir). 



(Al) 



The true scatter in equation lAll is very poorly known the- 
oretically, with values ra nging from d(logM vir) = 0.05 - 



0.1 d ex at fixed SFR (e.g., Button ct al. 2010; Finlato r et al 
2010 to 0.3 dex (e.g., iDe Lucia et al. 2006; Bcrton e et al 



20071 ). Our assumption of zero scatter is the simplest as- 



sumption possible. It allows us to study the evolution of 
galaxies at constant number density without reliance on 
theoretical expectations, and it allows the straightforward 
comparison to theoretical predictions for simulated galaxies 
selected in the same way. Moreover, this assumption is jus- 
tified by observations that find a tight relationship between 
UV luminosity and halo mass (e.g..|Gia valisco & Dickinson] 
I2OOII ; lAdelberger et al.ll2005l : iLee et al.ll2006l . 120091 ') . Even if 
there is appreciable scatter, we find evidence that the SFHs 
of galaxies may vary only by a scale factor (s ee jj 3), consis- 



tent w ith expectations from simulations fe.g.. lFinlator et al.l 
[2010). Therefore, even if our assumption of zero scatter be- 
tween UV luminosity and dark-matter halo mass is simplis- 
tic, the derived SFHs still should be reasonably accurate. 

Nevertheless, a non-zero scatter between UV luminos- 
ity and dark matter halo mass would affect the complete- 
ness in the descendants and progenitors of galaxies selected 
at constant number density at different redshifts. We study 
this effect using severa l models of galaxy evolution. For the 
iDe Lucia et al.l l|2006l ) SAM models we select galaxies at 
2 = 7.3 with SFR such that n(> 5') = 2 x 10"" Mpc"^ 
and track the descendants within the Millennium merger- 
tree at lower redshift directly. Because these models have the 
highest scatter at fixed SFR, d log Mvir = 0.3 dex, the com- 
pleteness in the descendants selected with constant number 
density at lower redshift declines. The bottom panel of fig- 
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ure IA2I shows that the completeness drops to 40% in these 
cases. However, we consider this model to be extreme as it 
has not been shown to reproduce the UV luminosity func - 
tions at hi gh redshift. In contrast, the iDutton et alj (|2010l ) 
SAMs and lFinlator et alJ (|2010l ) simulations have lower scat- 
ter, d(logMvir) ~ 0.1 dex, and these models broadly repro- 
duce trends in the SFR-mass and UV luminosity functions 
at high redshift. For these models we randomly adjust the 
mass of the dark matter halos in the Millennium simulation 
by the scatter and reevaluate the completeness. Figure IA2I 
shows that for these models the completeness in the descen- 
dants oi z — 7.3 galaxies remains >50% for all redshifts 
z ^ 3. In this case selection at constant number density se- 
lects the majority of the descendants as a function of redshift 
and has a negligible impact on our conclusions. 

Therefore, even if there is appreciable scatter in the 
SFR-halo mass relation, selection at constant number den- 
sity still tracks the descendants of the higher-redshift galax- 
ies, albeit with a higher contamination fraction. Further- 
more, we find that the majority of the the descendants that 
are "missed" by this selection as well as those galaxies that 
"contaminate" this selection (those without a progenitor se- 
lected at higher redshift) have SFRs within a factor of 3 
compared to the limit for constant number density. These 
galaxies appear to have evolutionary histories similar to the 
galaxies selected at constant number density and scatter into 
and out of the sample at each redshift. While this reduces 
the ability of the selection technique to track the direct pro- 
genitors and descendants of the galaxy population, it should 
introduce only a weak affect on the derived SFH. 

We therefore conclude that selection at constant num- 
ber density provides relatively fair, homogeneous sam- 
ples with which to study galaxy descendants and pro- 
genitors. The conclusion s here are similar to those of 
Ivan Dokkum et al.l (|2010f ) (see discussion at the beginning of 
this section). Comparing galaxy samples at constant number 
density seems preferable to study galaxy SFHs, compared to 
studies that use luminosity- or mass-selected galaxies. This 
is for the reason that the latter quantities evolve strongly 
in time. Samples selected at constant luminosity or stellar 
mass will be heterogeneous at different redshifts, potentially 
obfuscating any evolution in these quantities themselves. 



APPENDIX B: ON THE STELLAR MASSES 
DERIVED FROM MODELS V^ITH RISING SFRS 

One of the main conclusions from this study is that the 
SFRs of galaxies at constant number density increase from 
z = 8 to z = 3 and that this SFH is consistent with the evo- 
lution of stellar mass in these galaxies. However, their is a 
potential inconsistency with this later arguement as galaxy 
stellar masses used to make this comparison were derived 
from models using either co nstant or declining SFRs (see 
^ l2.3[l . lMaraston et al.l (|201Gr ) concluded recently that mod- 
els with rising SFRs provide better fits to the measured pho- 
tometry for z ~ 2 star-forming galaxies, and produce stellar 
masses larger than const ant or decaying r models by fac- 
tors of ~ 2 — 3 (see also iLee et al.ll2010bl ). Therefore, it is 
prudent to consider how the stellar masses of galaxies are 
affected using models with rising SFRs compared to those 
for constant or declining SFRs. 



We compiled multiwavelength photometry for galax- 
ies in the GOODS-S field with high fidelity sp < 
scopic ^ 



ectro- 



redshifts 3.0 < z < 6.0 (ICristiani e t al' '2000|: 
Vanzella et al.l l2008l : iPopesso et al.l 



Stanwav et al.l |200 



2009 : iBalestra et al.l |2010|) . The multiwavelength photom- 
etry was originally selected using the ACS zsso-band, and 
includes aperture matched colours in bandpasses from ACS 
BissVeoeiTTszsso, ISAAC JHK, and IRAC at 3.6, 4.5, 5.8, 
and 8.0 ji-m. The I SAAC and IRAC p hotometry are mea- 
sured using TFIT l|Laidler et al.ll2007f ) using the ACS zsso 
data as a prior for source positions. 

We derived stellar masses for this subsample of 
GOODS-S galaxies using decaying r m odels, ^ ~ 
exp(— i/r), following Papovich et al.l ( 2001 ). We used the 



2007 version of the iBruzual fc Charlod(|2003l 'l models with 
Chabrier IMF and fixed metallicity of 0.2 Zq. We allowed 
the stellar population age to vary from 10^ — 2 x 10^° yr, with 
r in the range 10^ — 10^" and r = oo (e quivalent to a constant 
SFR), with dust attenuation following ICalzetti et al.l (|2000f ) 
ranging from A{V) — 0.0 — 2.8 mag. We then derived stellar 
masses using the same models, but where the SFRs rise with 
time following equation |3] and where we fix the age (i) to be 
the lookback time from the galaxy redshift z to the fiducial 
formation redshift zf — 11. We allow for the same range of 
dust attenuation, and we allow the overall normalization of 
the model to vary in order to fit each galaxy because the 
evidence suggests that the SFHs of galaxies differ only by 
a fac tor that scales with mass (see § 3, and iFinlator et al.l 
I2OI0I ). We note that this SFH is the cosmological average 
of galaxies selected at constant number density. In reality 
individual galaxies will have unique SFHs that likely depart 
from this model. The point here is not to obtain necessarily 
the best representative fit to the photometry of each galaxy, 
but rather to understand how models with rising SFRs affect 
the derived stellar masses. 

Figure [Bll shows several examples of galaxies with stel- 
lar masses that would fall in a sample selected at constant 
number density, n — 2 x 10^* Mpc^''. The curves in each 
panel show the rising SFH and declining SFH models with 
the "best-fit" parameter values from the model with the 
minimum reduced x^ value with parameters listed. In all 
cases the stellar population ages and stellar masses have 
larger values for models with rising SFHs. The minimum re- 
duced x^ are close in all cases. This is a result of the fact that 
the SFH is poorly constrained by modeling the broa d-band 
photometry of galaxies (e.g., IPapovich et al.ll200ll ). How- 
ever, in general the models with rising SFHs provide fits 
that are statistically indistinguishable as models with de- 
clining models, which is noteworthy given that here we have 
tested only a single rising SFH (with 'l'(i) ~ t^'""^) compared 
to the declining models, which spanned the full range of e- 
folding time-scales (see above). Finally, we remind the reader 
that we expect models with rising SFHs to represent average 
galaxies with ongoing star-formation. Realistically, galaxies 
have individual, stochastic SFHs. In particular, at lower red- 
shifts, 2 < 3, galaxies will exhibit "suppressed" SFRs, and 
these are consistent with declining SFHs. The rising SFHs 
advocated here corresponds to the average evolution of the 
star- forming galaxy population at early times (3 < « < 8). 

The results of our fitting analysis show that stellar 
masses derived with rising SFRs generally produce higher 
stellar masses than decaying r models, and that the increase 
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Figure Bl. Fits to tlie broad-band piiotometry for several galax- 
ies from the spectroscopic GOODS-S subsample. Each panel 
shows the results for an individual galaxy. The ACS, ISAAC, and 
IRAC photometry are indicated by yellow squares. The curves 
show the best-fit models for models with rising SFHs (brownish- 
red curves) and declining SFHs (light-blue curves) with best-fit 
parameters inset in each panel. The full range of model param- 
eters used are listed in the text. In general, models with rising 
SFHs yield higher stellar masses and older ages, but both models 
have best-fit reduced y^ that are statistically generally indistin- 
guishable. 
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Figure B2. Stellar masses derived from decaying r models where 
Vt ~ exp(— t/r) versus those derived using the rising SFH derived 
here (equation jsj . The panels show galaxies with spectroscopic 
redshifts from GOODS-S in two redshift bins of 3 < z < 4.5 and 
4.5 < z < 6. Representative error bars are shown. The black, 
short-dashed lines show the unity relations. Stellar masses de- 
rived with rising SFHs generally produce higher stellar masses 
than decaying r models. The increase in stellar mass generally 
increases with time. The stellar masses derived from models with 
rising SFRs for galaxies at 3 < z < 4.5 are a factor of 1.6 higher 
(indicated by the red, long-dashed line, bottom panel), while the 
offset is smaller (factor of 1.1, red long-dashed line, top panel) for 
galaxies at 4.5 < z < 6.0. 



ies at 2 ~ 3 — 4 may need to be revised upwards by fac- 
tor of order 2, simila r to the conclusio ns of JMaraston et al.l 
(J2OIOI ) (although see JLee et ai]l2010bl . who find that rising 
SFHs overpredict stellar masses of galaxies in simulations 
by factors of up to ~2). The masses of galaxies at z > 4 are 
mostly unchanged. Moreover, this increase in stellar mass 
from rising star-formation histories for galaxies at z < 4 is 
approximately the amount needed to reconcile the offset be- 
tween the measured SFRs and stellar masses (figure [3]). If 
the stellar masses of the UV-selected galaxies are higher by 
0.2-0.3 dex, it would adjust them closer to the derived SFH. 

This paper has been typeset from a T^^jK/ ffl^^X file prepared 
by the author. 



in the stellar masses generally increases with the passage of 
time. Figure [B2l comDares the stellar masses derived from the 
decaying r models to those derived using the rising SFRs, 
divided into two subsamples of redshift, 3.0 < z < 4.5 and 
4.5 < z < 6.0. The stellar masses derived from rising SFHs 
for galaxies at 3 < 2 < 4.5 are a factor of 1.6 higher on aver- 
age, while the offset is smaller (a factor of 1.1) for galaxies 
at 4.5 < 2 < 6.0. Therefore, if the model SFH posited in 
§ 3 is correct, then the stellar masses of star-forming galax- 
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